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The stability of the shell of wrapped D6-branes on K3 
is investigated from the point of view of supergravity. We 
first construct an effective energy-momentum tensor for the 
shell under the reasonable conditions and show that super- 
symmetric solutions satisfy Israel's junction conditions at ar- 
bitrary radius of the shell. Next we study the perturbation 
of the whole system including the self-gravity of the shell. It 
is found that in spite of the existence of wrapped D6-branes 
with negative tension, there is no eigenmode whose frequen- 
cies of the shell and the fields are imaginary numbers, at any 
radius of the shell. Furthermore, when the radius of the shell 
is less than the enhangon radius, resonances are produced, 
and this indicates a kind of "instability" of the system. This 
can even classically explain why the shell is constructed at 
the enhangon radius. 



I. INTRODUCTION 

One of the central issues in quantum gravity is the 
resolution of singularities, i.e., the initial singularity in 
our universe or a central singularity in black holes. It 
is expected that string theory, which is one of the most 
prominent theories including the quantum gravity, may 
resolve the issue of such singularities. From the gen- 
eral relativity theory point of view, the appearance of a 
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naked singularity is a serious problem because we can- 
not predict what happens in the future. Therefore, it is 
particularly important to study whether the naked singu- 
larities are excised in string theory or not. An interesting 
example of the naked singularity is repulson singularity, 
which appears in supersymmetric solutions with spheri- 
cal symmetry in supergravity theory [0-0] . Kallosh and 
Linde showed that all test particles, either massless or 
massive, cannot touch the repulson singularity. This par- 
ticular type of singularity is a globally naked singularity 
and the weak cosmic censorship [|J seems to be violated. 
This is essentially different from the extreme Reissner- 
Nordstrom black hole where a locally naked singularity 
exists inside the event horizon. 

Recently, Johnson, Peet and Polchinski analyzed the 
motion of a D6-brane probe wrapped on a K3 manifold 
in the type IIA supergravity solution and proposed the 
enhangon mechanism which excises the repulson singu- 
larity HH . The higher curvature corrections on the world 
volume of D6-brane induce a negative D2-charge and a 
negative D2-tension (see Refs. J7|— p~0[ for details), and the 
total tension of the wrapped D6-brane probe vanishes at 
a special radius r — r e larger than the radius of the re- 
pulson singularity |ra|. This fact led the proposal that 
gauge symmetry is enhanced at the radius r — r e , which 
is called the enhangon radius, and that a shell consisted 
of the wrapped D6-branes is located at the enhangon ra- 
dius. As a result, the original repulson singularity can be 
removed due to the flat geometry for r < r e . 

Since there is no precise dual theory to this type of 
supergravity theory at the present stage, this motivates 
us to investigate the picture of the enhangon mechanism 
from supergravity side more deeply. It is especially inter- 
esting to examine whether or not the enhangon geome- 
try is a solution of the supergravity theory including the 
self-gravity of the branes, and whether the shell of the 
wrapped D6-branes is stable or collapses by the gravita- 
tional effect. 

So, in this paper, we first construct an effective energy- 
momentum tensor of the shell consisted of the wrapped 
D6-branes and study the static solutions by using Israel's 



junction condition ]ll| . Next we study the spherically 
symmetric perturbation of the whole system and ana- 
lyze the stability of the obtained solution. By numerical 
analysis, it is found that there is no eigenmode whose fre- 
quency lo of the shell and the fields is imaginary. To our 
surprise, although the tension of the wrapped D6-brane 
is negative in the region r < r e , the shell of the branes 
never has eigenmodes with u> 2 < if we take into account 
the interactions between the shell and the fields, which 
include the self-gravity of the shell. 

Furthermore, it is shown that when the radius of the 
shell is less than r e , resonances are caused at proper 
oscillation frequencies of the shell, namely the system 
becomes "unstable" in this region. On the other hand, 
when the radius of the shell is greater than or equal to r e , 
the system is really stable, that is, there is no eigenmode 
with u! 2 < and no resonance is produced. These facts 
indicate that it may be able to explain even classically 
the dynamical reason why the shell is constructed at the 
cnhangon radius. 

The outline of this paper is as follows. In Sec. y, we 
briefly review the D6-D2 brane solution and show that 
the velocity of the D6-brane probe reaches the speed of 
light at the enhangon radius. In Sec. Ill, we derive the 
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effective energy-momentum tensor of the shell of the D6- 
D2 branes under the reasonable conditions, and give the 
energy conservation law and the dilaton equation on the 
Gaussian normal coordinates. In Sec. IV, we show "the 



enhangon geometry" with the shell at any radius is the 
solution of the supergravity based on the Israel's junc- 
tion condition. In Sec. [v|, we first derive perturbed equa- 
tions of the gravitational and matter fields inside the shell 
located at arbitrary radius. Next, we derive perturbed 
junction conditions of the fields on the shell, the equation 
of motion of the shell and perturbed equations of fields 
outside the shell. Finally, using numerical analysis, we 
study the stability of the shell and the fields. Sec. [Vl] is 
devoted to conclusion and discussion. Throughout this 
paper, we set a' = 1. 



II. ENHANgON GEOMETRY AND THE 
RELATIVISTIC MOTION OF A D6-BRANE 

The low-energy IIA supergravity action is given by 
1 
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dC( 3 ) +C(i) A dB( 2) , where 
C( 7 ) is the dual field of Cm. The hats on the fields and 
the differential denote that the contractions are done with 
the string metric Gmn- 

Let us consider a D6-D2 system in which there are N 
D2-branes and N D6-brancs. The supersymmetric solu- 
tion of the system is 

ds 2 = G M Ndx M dx N 

= ZpZ^rjpv dx»dx v + z\zl& l3 dx l dx j 
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C (7) = V(gZ 6 )- l dx° A dx 4 A dx 5 A dV {K3) , (2.5) 

where /i, v = (0,4,5), i, j — (1,2,3), and p, q — 

(6, 7, 8, 9). The G pq denotes the metric of the K3 mani- 
fold with unit volume^. The D2-branes are aligned along 
the 4, 5 directions, while the D6-branes are aligned to the 
4, 5, 6, 7, 8, 9 directions. 

The harmonic functions are expressed as 
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dVtKS) is the volume form of the K3 manifold of unit vol- 
ume. To keep the tension of a wrapped D6-brane, which 
will be investigated in detail later, positive at infinity, we 
assume that V > V* := (2-7r) 4 . 

Because of the negative charge of the D2-branes, there 
is a repulson singularity at r = \r 2 \, where Ricci scalar 



curvature diverges as (r — \r%\ 



-3/4 



To investigate the 



causal structure of this geometry with the repulson sin- 
gularity, let us solve a radial null geodesic in this solution. 



By using Eq. (2.2), the future- ingoing null geodesic be- 



haves near the singularity as 
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where t := x°. This means that t is finite when r — > \r 2 \ 
and hence the repulson singularity is a timelike singular- 
ity, as depicted in Fig. 1. To resolve the repulson sin- 
gularity, Johnson, Peet, and Polchinski have proposed 



where Gmn is the metric in the string frame, $ is the 
dilaton field. H^ is the field strength of the two-rank 
anti-symmetric tensor field B( 2 ) ■ F(p+2) is expressed by 
the (p+l)-form Ramond-Ramond (R-R) field C^ p+ i) as 



' Although we cannot write down the metric explicitly, it 
is well known that the Ricci curvature is zero. For later 
convenience, here, we shall formally denote the metric as 
G^ 3) dx v dx q 



the enhangon geometry which is flat inside the unique 
enhangon radius |5j]. 
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Repulson Singularity 



FIG. 1. The Penrose diagram of a D6-D2 supergravity so- 
lution is shown. The dashed line represents an orbit of a 
D6-brane wrapped on the K3, while the dotted line repre- 
sents the orbit of an ingoing-radial null geodesic. 

If we ignore higher order corrections including the cur- 
vature corrections, the action of a Dp-brane is given by 
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where fi p = (2n) p is the charge of the Dp-brane and g a t, 
is the pull-back of the string metric Gmn'- 
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B a t, is the component of the spacetime NS-NS fields S/ 2 ) 
parallel to the Dp-brane and F a b (and F) is the field 
strength of the gauge field living on the Dp-brane. In 
this paper, we consider the D6-D2 brane system with 
vanishing B( 2 ) an d Fob, so that the action becomes 



jDp 



-lip dP +i ^e-' 

+M P I dP +1 tC (p+1) . 



det g ab ) 



(2.11) 



Now, consider the motion of a wrapped D6-brane on 
the K3 manifold. If we take the higher-curvature cor- 
rection terms in the K3 to the action of the D6-brane 
into account, a D2-brane with a negative D2-charge is 
induced (see Jz| p~0|| for details). Hence, the total action 
of the wrapped D6-brane on the K3 becomes S® 6 — S® 2 . 
After integration of the wrapped coordinates on the K3, 
we obtain the effective action of the brane which is just 
a membrane: 
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where V(r) is defined by V(r) := VZ 2 /Z 6 . Note that the 
tension of the brane vanishes at the enhangon radius r e 
satisfying the following relation: 
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The explicit form of the enhangon radius is 
2V 



V - V. 

If we take a static gaugepl 
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the reduced Lagrangian density L is 



* Since D2 brane is a compact surface, the space of x 4 and x 5 
should be also compact in the static gauge, to be exact. Here, 
we take this space as two dimensional torus of unit volume. 
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where dots denote the differentials with respect to t. 

The energy conservation law gives us the following 
equation 
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where E is a constant. If the probe has an ingoing- 
velocity at r — ri 3> r e , E is strictly positive. If the 
brane probe stops at a radius r m > r e , then the left 
hand side would be zero at r = r m . This leads to con- 
tradiction because E is positive. Hence the brane probe 
cannot stop outside of the enhangon radius r e . At the en- 
hangon radius, the first term in the square bracket should 
diverge infinitely because V/Zq — V*/Z2 — at r = r e . 
This means that the velocity of the brane probe reaches 
the speed of light at the enhangon radius (See Fig. ||)! 

It seems strange at first sight that according to the 
picture by Johnson, Peet and Polchinski in Ref. g, the 
branes should collect near the enhangon radius and that 
they cannot leave the location of the enhangon radius, 
while we have just seen, the probe discussed above moves 
at a speed of light near the enhangon radius. Hence we 
need further analysis about this enhangon radius and the 
mechanism beyond the probe approximation. In next 
sections, we study the whole system including the self- 
gravity of the branes and the self-interactions between 
the fields and the branes. 



in the neighborhood of the D6-D2 branes shell such that 
it moves along \ — Xo constant. 

Here after, we shall impose spherical symmetry for sim- 
plicity. Hence, we can choose the Gaussian normal coor- 
dinates (x,C a ) as X — x( r it) an d C° = C°( r )*)> while the 
other coordinates are same as 6, </>, x m (m = 4,5), x p . 
The nine-dimensional embedded metric is rewritten by 

^ 2 | x =const. =<la b dCd( b 

= -qoo(d( Q ) 2 + GeedO 2 + G H d$ 

+G mn dx m dx n + G pq dx p dx q . (3.3) 

To derive the effective energy-momentum tensor com- 
ing from the wrapped D6-branes, we shall impose the 
following conditions: 

1. we can take a static gauge such as £° = £°, £ m = 
x m and £ p = x p for each D6 brane and £° = (° and 
£m _ x m £ or eacn D2 brane. 

2. each D6 (D2) brane moves along the radial direc- 
tion only. 

3. each D6 (D2) brane is distributed uniformly on the 
shell. 

The condition 2 implies, in other words, that each D6-D2 
brane moves on the x ~ Xo: ( a — const (a^O). 

In order to see the way how to derive the effective 
energy-momentum tensor of the D6-D2 brane shell, it 
would be enough to discuss J2 i=1 SP? only. Using the 
conditions 1 and 2, the action can be reduced as follows: 



III. THE ENERGY-MOMENTUM TENSOR OF 
THE D6-D2 BRANES SHELL 

By the conformal transformation such as Gmn — 
e _ */ 2 GM"jv, we obtain the low-energy effective action in 
the Einstein frame: 
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where Gmn is the Einstein metric and all the contrac- 
tions are done with Gmn- i represents the number- 
ing of each wrapped D6-brane of which a shell consists. 
In p^Jl^ l , the energy- momentum tensor is obtained when 
the shell is fixed at a radius. In this section we will ob- 
tain the tensor by the following covariant way so that it 
can be applied to more general cases. 

In general relativity, it is well known that there is the 
following local (Gaussian normal) coordinates 



ds 2 = d X 2 + q ab dCd( b 
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Because of the condition 3, each D6-brane is distributed 
uniformly over the two-sphere (9,(f>). Introducing the 
number density p e satisfying 



dOd^iGeeG^)* p e = N, 



(3.5) 



we can rewrite Eq. (3.4) as the following ten-dimensional 
form: 

S 6 = -Me / d 9 Cei' s '(-dctq ab )ip 6 

= - M6 / dx w S( x - Xo ) e l*(- det G)ip 6 . (3.6) 



Note that pe is a function of (° only. As concerns the 
subspace (x m ,x p ), we can easily calculate the energy- 
momentum tensor by using the usual formula 
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where C m is a matter Lagrangian. In the subs pace 
(C°, 0, <p), however, we should imp ose the constraint ( |3.5| ) 
when we derive T ab by using Eq. (|3.7| ). Let us define the 
three-dimensional unit timelike flow vector along each 
D6-brane by V a = (V°, V e , V*) := {y/-q oa , 0, 0). Then, 
in terms of the current vector J a = p e V a , Eq. ( |3.5[ ) is 
rewritten in the following covariant form: 
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where ; denotes covariant derivative with respect to the 
subspace x a — ((°,9,(f>). So, as seen on page 70 of 
Ref. fPl, 



Sp 6 = ^-(V a V + q a0 )Sq a ^. 

Therefore, combining with 5y/—q/8q a P 
we get the tensor components as 
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(3.10) 

Just like the case of D6 branes, ^j=i ^H) can be a ^ so 
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where pi is a number density over the six-dimensional 
space, K3 x S 2 . It satisfies a conservation law as 



(3.13) 

It is interesting to note that the energy-momentum 
tensor has several particular features. Firstly, the dom- 
inant energy condition is violated near the enhangon 
radius r = r e because X^ 2 + T^ 6 is proportional to 
(V/Z 6 - V*/Z 2 ) ~ and hence \T§? + T Q f\ « \TP q 2 + 
T^ 6 \, where p, q implies an orthogonal basis. Sec- 
ondly the energy-momentum tensor looks like dust in 
the three-dimensional subspace (x ,9,(f>), while it looks 
like a varying cosmological constant in the other subspace 
(x°,x 4 ,a: 5 ), due to along which directions D6-D2 branes 
are aligned. 



IV. ISRAEL'S JUNCTION CONDITION AND 
ENAHNgON GEOMETRY 

If the shell sits on a radius r = tq (> r e ), timelike 
hypersurface x = Xo corresponds to r = rn one. On the 
hypersurface, Israel's junction equations |11| are 
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where S D b 6 and S D2 are functions in front of (5-function 
of T^ b 6 and T^ 2 , respectively. Using the unit normal 
vector % = {dx)Mi the second fundamental form K ab 
is defined by K ab :— n a -b- [f]± simply denotes /+ — 
/_, where / + and /_ are quantities evaluated on the 
outside and inside of the r = r$ timelike hypersurface, 
respectively. 

According to the enhangon mechanism, the inside of 
the shell is replaced with a flat geometry. Reminding 
that Gmn = e~ s'Gmtv, we obtain 



dx = Z 2 1B Z 6 16 dr. 
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Thus, we can easily calculate all components of K ab at 

any r = tq : 
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dPxdOdcj) (GggG^ det G pq )i p 2 = N. (3.12) 

Defining the seven-dimensional timelike vector U" = 
(U°,U e ,U<i',U p ) := (y/-q oo ,0,---,0), the effective 
energy-momentum tensor is 
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We can check that Eq. (4.1) is s atisfied at a ny radius 
r = ro with the help of Eqs. ( |3.5D and (3.12). So, we 
obtain a whole static solution without any repulson sin- 
gularity, taking the self-gravity of the shell into account. 
This is one of the main results in this paper Jig ] . 

Next, let us consider the whole energy-momentum ten- 
sor Tmn coming from both supergravity and the shell 
parts and derive the energy conservation on the shell. In 
terms of the unit timelike vector u M on the shell, the en- 
ergy conservation is written by u m W n Tmn — 0. To de- 
rive the conservation equation effectively on the shell, we 
shall put a step function 0(x~ Xo) m front of the matter 
Lagrangian density of the supergravity part in Eq. ( |3.1| )p|, 
i.e., / d 10 xV^GC m -► / d 10 xV^Gd(x-Xo)£ m , where 
C m = -(V$) 2 /2-e*/ 2 F 4 2 /2-4!-e- 3 */ 2 F 8 2 /2-8!. Then, 
Tmn is derived by the action ( |B.l| ) as follows: 
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f(x) means that f(x) is integrable and J °_ f(x)dx 



0(e). Integrating Eq. ( 4T4 ) from x = Xo~ e to x — Xo+e, 
we obtain the following equations: 
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As easily seen, this equation is always satisfied when the 
shell is static. 

Finally, let u s consider the dilaton equation. Varying 
the action (3.1) with respect to <&, we get the dilaton field 
equation 
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In terms of Gaus sian normal coordinates ( |3.2j ), the left 
hand side of Eq. (4.16) is expressed as 



V 2 $ = 
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Since th e seco nd term is a regular function of x, integrat- 
ing Eq. ( |4.16| ) from x = Xo - e to x = Xo + e, we get a 
junction condition for the dilaton field such as 
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In the following sections, we will consider the pertur- 
bation of the solution. 
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and S m ~ D2 MN ■= S D6 mn + S D2 MN - Note that 
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step function, all components of S^mn and S dllaton MN 



are smooth functions of x an d they can be calculated on 
the outside of the shell. Then, u ai W n Tmn = is written 



by the Gaussian normal coordinates ( |3.2j ) as 
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where D denotes covariant derivative with respect to the 
metric qMN = Gmn — ^m^at- Here, a regular function 



^We thank to Kouji Nakamura for suggesting this. 



V. STABILITY ANALYSIS OF THE D6-D2 
BRANE SHELL 

To avoid complication, let us focus on the following 
perturbed metric in Einstein frame. 



ds 2 = 
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Z 6± (t,r) = Z 6± (r)[l + SZ 6± (t,r)], 



(5.2) 



where ± denotes outside and inside of the shell. / denotes 
an unperturbed background field /. The dilaton field 
$± and Ramond-Ramond fields C@)± and C(j)± are also 
written by 



C 3 ±(t,r) = C 3 ±(r)+5C 3 ±(t,r), 
C 7± (t,r) = C 7± (r)+6C 7± (t,r). 



Eq. (5.9) indicates 



(5.3) 



Here the fields C 3 and C 7 are components of the Ramond- 
Ramond fields CW) and C(j\ respectively: 



C( 3 ) = C 3 dtf\dx 4 f\dx b . 

C (7) = VC 7 dt A dx 4 A dx 5 A dV(#3). 

The above unperturbed fields are defined as 



Mr) 
Z 6 {r) 



Z_2-{r) = 1 + r 2 /r 
Z 2+ (r) = l + r 2 /r 

Z 6 -{r) = 1 + r 6 /r 
Ze+(r) = l + r 6 /r 



(5.4) 
(5.5) 



(5.6) 
(5.7) 



and <&±(r), C 3 ±(r) and C 7 ±(r) are given by Z 2 ±(r) and 
Z 6 ±(r). 

Assume that the shell sits on r — ro before it is 
perturbed. Since the inside of the geometry is just 
flat, all functions /_ are constants determined by the 
continuity of the outside and inside of the metric, i.e., 
/-=/+(r = ro). 



A. The perturbation of the flat geometry inside the 
metric 

The perturbed energy-momentum tensor STmn— does 
not include any linear terms of 5$-, SC 3 - or 5C 7 - be- 
cause the background energy-momentum tensor Tmn- 
is zero. Hence, the perturbed Einstein equation is sim- 
ply SRmn- = 0. Defining derivative operator V 2 _ :— 
9(9* + did 1 , the equations are written by 



V 2 L{6^-+5Z 2 -+6Z 6 -) 

-df (8 <5tf_ - 6 5Z 2 _ + 2 5Z 6 _) 

d t d l {2S^- - 5Z 2 _ + SZ 6 -) = 0, 

V 2 . (<5#_ + <5Z 2 _ + bZ§J) = 0, 

Vi (m- - 8Z 2 _ - 5Z 6 J) Sn 

+d l d j (8 £*_ - 2 5Z 2 _ + 6 <5Z 6 _) 



0, (5.8) 



(5.9) 



(5.10) 



0, (5.11) 



If we take 



SZ 2 



SZ K 



0. 



S*± = $±(r)e iut , 

5Z 2± =( 2± (r)e^, 

SZ 6± = ( 6± (r)e i « t , 

S$±=r ] ±(r)e i " t , 



(5.12) 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



e- = 2(C2--c 6 - 



(5.17) 



Substituting this into Eqs. (5.8)-(5.12), we obtain the 
following equations 



V*<?*_ = 0, 



£- = (2- 



<e- 



(5.18) 



(5.19) 



The pertu rbed field equation of the dilaton is same as 
Eq. ( paj ); 



Vi<5$_ = 0. 



(5.20) 



Imposing a re gular boun dary condition at r = 0, the 
solution of Eqs. ( 5.18 ) and ( 5.20| ) are 



Ur)=A, 



r/_(r) = A 



sin(Jlr) 
fir 



sin(Jlr) 
fir 



(5.21) 



(5.22) 



where 51 = \] Z 2 -Z§- u>, A<y and A§, are integral con- 
stants. Note that when the shell is absent, the frequency 
u) must be real. Otherwise, i.e., when u is imaginary, 
the eigenfunctions diverge at infinity. However, when the 
shell exists at some radius, the gradient of eigenfunctions 
becomes discontinuous at the shell and there is a possi- 
bility that the exponential growth of the eigenfunction 
may be changed by the exponentially decaying function 
at the shell. In this sense, the stability of the shell is not 
trivial. 



B. The perturbation of the junction conditions on 
the shell 

Suppose that th e orbit of the shell is r = r^ + SR(t) 
in the metric (5.1). In this and next subsections, we 
mainly discuss the values around the shell. Hence all the 
variables are evaluated on the shell if we do not mention 
in these subsections. We also drop the index ± of the 
continuous variables on the shell. Then, the (linear or- 
der) continuity of the metric and dilaton field yields the 
following equations, 



<5$+ 


= £$_ 


- $V 8R, 


SZ 2+ 


= 5Z 2 _ 


^2 


5Z 6+ 


= 5Z 6 - 


-^SR, 



(5.23) 



where a prime denotes differential with respect to 



Firstly, let us solve Eq. (4.15) and show that the dila- 
ton field $ is constant along the shell, i.e., d$/d(° = 0. 
The components of the vectors u M = u M 
Sn M are expressed as 



8u M and 



n M _ n M 



n M = (n*,n r ,0,---,0) 

= (o,^*z 6 -*,o,---,o) 

Sn M = (zfzpSR, 

~Z2*Z^(8*-6Z 2 -6Z 6 ),0,- 



,0) 



^M 



= (z/°z 6 i6 ,o,o,---,o), 



Su M = i^-zi 6 z^ {m + sz 2 + sz 6 ), 



Z^Zi°5R,0,---,0 



(5.24) 



where dot means derivative with respect to t. After 
straightforward but long and tedious calculation, we get 

<y(£ daaton „„u JI 
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(5.25) 



with help of Eqs. (5.19) and (p.24). Finally, substituting 



Eqs. ( |5.25| ) and fl5.29|)into Eq7([llJ), we get the follow- 
ing simple equation 



<5$_ = 0, 



(5.30) 



which means that the dilaton field is constant along the 
shell in the first order of the perturbation. It should 
be noted, however, that it does not necessarily imply 
6$>-(r) — 0. For example, the ge neral solution of the 
perturbed dilaton field equation is ( 5.22). W hen the fre- 
quency O (or uj) is real, the condition ( 5.3C ) means that 
the <5$_ (r) becomes node at r — tq. The frequency of the 
perturbation becomes discrete to satisfy this condition 
when the perturbation of the dilaton field has nontrivial 
configuration. Of course, there may exist a solution with 
A$ = 0. In this case, the frequency shows continuous 
spectrum. On the other hand, when m is imaginary, the 
eigenfunction grows exponentially with r, and it cannot 
become zero unless ^4$ = 0. Hence <5<J>_ = for the 
imaginary frequency, i.e., unstable modes. 

Now, let us focus on the perturbation of the junction 
equations (4.1). It is noteworthy that the right hand 
side does not include any 8<$>- terms, pr ovid ed that we 
calculate it by using the inside metric (5.1), as shown 
above. So, all non-vanishing equations are: 
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zl 



3T2 

Z 2 
SR, 



Z 6 



SZ, 



(5.31) 



r/cRR M N\\ n 

d(S mnu n )L =Yn = 0, 



(5.26) 



where we used Eq. (5.23 ) to d erive the second equation. 
The first term in Eq. ( 4.15 ) is rewritten by 
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-SZ' a _ + ^%SR, (5.32) 



8r, 



3^ + ^)^.(5.33) 



By solving these equations with respect to S^>' + , SZ 2+ , 
and SZ' &+1 we obtain 
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(5.34) 



This 



W = pee' 
quantity 



y*pi e" 



(5.28) 



can be expanded as 

+ 5{u M \7 N S m - D2 MN ) in terms 



u M V N S D 6 -D2 MN \ r=ro 

of both the metrics outside and inside the shell. There 



fore, it is convenient to expand Eq. (5.27) in terms of the 
metric inside the shell. As a result, we get 
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(5.35) 



SR. (5.36) 



The left hand side in Eqs. (5.34)-(5.36) are interpreted 
as a difference between inside and outside quantities. So, 
if the number densities p 2 and p% of D2 and D6 branes 
on the shell approach zero, the difference also approaches 
zero. 

Similarly, we can obtain 5< &j_ fr om the junction condi- 
tion of the dilaton field Eq. (4.1 8|) as follows: 
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We also have the following equations from the RR field 
equations using Eq. (|5.23): 
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SC' 3+ 
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l23 SR, (5.38) 
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2gr 2 Z 2 



35$_ + 5$_ 
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4 ^SR. (5.39) 



3 r o z 2 



C. Equation of motion of the shell 

When the locus of the shell made of the wrapped D6- 
branes shifts from r — r$ to r = R(t) = ro+5R(t), where 
|5i?(£)| <C tq, the effective action of the shell up to the 
second order of SR(t) is given by 
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(5.40) 



From the action ( 5.4C ), we have the following equation 
of motion of the shell: 



p 6 VZ 2 (R)-p 2 Z 6 (R) 



5R 



4Z 6 (R) 



z' e (R) 

Z e (R) 



(35$ - 75$ + SZ 2 - 7SZ 6 
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/i 2 



4Z 2 (R) 



Z 2 (R) 



Z 2 (R) 

+ (<5$ + 35$ + 3SZ 2 + 3SZ 6 y 
+gfi 6 V5C! r - g» 2 5C 3 , 



(5$ + 35$ + 35Z 2 + 35Z 6 ) 



(5.41) 



where the primes denote differentials with respect to R. 
8C 3 and SCj are the perturbed parts of the Ramond- 
Ramond fields respectively. 



The fields in Eq. (5.41) are discontinuous at the shell 
r = R(t\ so that the equation needs to be averaged near 
the shell. Recalling that the average of A(R) at the shell, 
{A}, is defined by {A} = h(A- + A+), the equation of 
motion of the shell up to the first order of the perturba- 
tion can be written as 
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+ |m6^(5C 7 _ + SC 7+ y - |a* 2 (5C 3 - + SO*)', 






(5.42) 



where we can take Z!, 



-r 2 /r 



2 yl 
Z 6 



~ r 6/ r o 2 - Here 
we have also used Z^_(R) = Z! 2 _(R) = 0. Using the 
relations between the inside and the o utside quan tities 
near the shell, that is Eqs. (|5 19[ ) and ( |5.34| )-( |5~39| ), we 
finally obtain the equation of motion of the shell: 

ASR = B + C6R, (5.43) 

where the definitions of A, B, C are as follows 

A = p 6 VZ 2 - p 2 Z 6 (5.44) 
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(5.45) 
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Here, we have used SC 3 _ = SC 7 _ — inside the shell. 
Moreover, using a := V/V* and x := ro/|r2|, A, B and C 
have the following simple expressions: 

x 



A= — [{a - I )x - 2a] 



(5.47) 
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A*2 / 7a 



128r 2 x V x + a x — 1 
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(5.48) 
(5.49) 



D. The perturbation of the geometry outside the 
shell 

The calculation of the perturbation of the outer region 
is tedious but straightforward. First we examine the mat- 
ter field equations. The perturbed equations of the RR 
fields become 

SC 3+ ' = -^6Z 2+ - <5*+ + 5$+) C' 3 , (5.55) 

SC 7 + = -^{4SZ 6+ + 3S^+ - 3<5$ + ) C' 7 . (5.56) 



5R = 6Re tuJt (5.50) 

B = Be tu> \ (5.51) 

the equation of motion of the shell becomes 

aSR = -B (5.52) 

a = uj 2 A + C. (5.53) 

Therefore, we obtain the amplitude of the forced oscilla- 
tion: 



6R 



B 

a 



(5.54) 



Using the above equation, we can investigate whether 
resonances are produced by the forced oscillation or not. 
If a can be taken 0, then the amplitude SR becomes in- 
finity and a resonance is caused. This indicates a kind 
of "instability" of the shell within our perturbation the- 
ory. F urthe rmore, from the numerical analysis of the 
Subsec. VE, the square of the frequency io cannot take 
negative values. This fact and the behavior of the a, tak- 
ing notice that C > at any radius, lead us to classify 
the region of the radius of the shell into the following 
cases. 

(i) When r e < r , A > (positive tension) and C > 0, 
so that the resonance is not caused, 
(ii) When r = r e , A = (tensionless) and a = C > 0, 
so that the resonance is not caused. 
(iii) When r < r e , A < (negative tension) and C > 0, 
so that there is necessarily a real frequency w satisfying 
a = at an arbitrary radius of the shell. Therefore, a 
resonance is produced and the shell oscillates with huge 
amplitude at this radius. In that sense, the shell is "un- 
stable" inside the enhangon radius. 

From the above analysis, the "unstable" region of the 
shell, where resonances are produced, is r^ < r e . If the 
shell of the wrapped D6-branes enters the inside of the 
cnhanqon radius r e , the shell becomes "unstable" because 
of the occurrence of the resonances under small fluctu- 
ations of the fields, and the shell may push back to the 
enhancon radius, where it is stable. This can explain why 
the shell is constructed at the enhangon radius. 



The R-R fields are expressed by the perturbations of 
other fields and obtained just by integrating these equa- 
tions after we determine other field. The perturbed equa- 
tion of the dilaton field becomes 
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Next we turn to the gravitational equations. The first 
order perturbation of the Ricci curvature becomes 
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On the other hand, the energy momentum tensor is 
perturbed as 
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Here we have used Eqs. ( ^5.55 ) and (5.56). 

Since STqq + + ST mm+ = 0, combining the (00) and 
(mm) components of the Einstein equations gives the 
following useful relation 
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By assuming the time dependence of these pertu rbed 
variables as harmonic oscillator as Eqs. ( 5.15 )-( 5.16 ), we 
find 
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By the (0,0), (r,r) and (p, q) components of the per- 
turbed Einstein equations, we find the equations of £2+ 
and C6+ 1 
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By using the relation ( 5.69| ), the dilaton field equation 
becomes 



V = *+&+' + &+') 



% (ft + ft ' --" 



3/ 7/2 7/2 
/ ^2+ ^6+ 



4 V Z 2 + ^6+ 



(C2+ - C6+) - W 2 ^2+^6+»/+- 



(5.72) 



E. Numerical analysis 

Now, everything has been prepared for the pertur- 
bation analysis, i.e., the pert urbed solu tions of the ge- 
ometry inside the shell ( 5.19) an d ( |5.2l|), th e mat ching 
and the junction conditions ( 5.23) an d ( |5.34 )-( |5.37 ), the 
equation of motion for the shelT(5.43), and the perturbed 
equations in the outer geometry (5.69)-( p.72j ). Since we 
know the solutions inside the shell, we can integrate from 
the To to outside region with the boundary value and 
the junction conditions of the perturbed functions at the 
shell. If there is an eigenmode with a negative eigenvalue 
w 2 , the background solution is unstable against the per- 
turbations. 

Our numerical calculations show, surprisingly, that 
no eigenmodes with negative u 2 can be found for any 
radii (ro > \r 2 \) of the shell. Furthermore, as shown in 



Subsec. |VC| , we found resonances for the case (hi). The 
followings are two typical examples for the numerical cal- 
culations. 

Fig. pi shows configurations of the stable modes when 
the shell is located at r = rn > r e . The enhancpn radius 
of the parameters in Fig. g is r e = 10.8123. Outside 
of the shell, the perturbed functions couple each other 
complicatedly and they show disordered behaviors. The 
stable modes show a continuous spectrum as usual. 

When the shell is located inside of the enhangon radius, 
the configurations of the perturbed functions are similar 
to the previous case r > r e except for the existence of 
the resonance radii. For the same parameter in Fig. 0, 
the resonance radii are r reso — 2.74099, 10.8071 in \r 2 \ < 
r reso < r e , where \r 2 \ = 2.59758 and r e = 10.8123. Fig. | 



shows configurations when the shell is located near the 
resonance radius. The perturbed functions becomes very 
large at the shell by the resonance. 

When the shell crosses the enhangon radius, the total 
tension of the shell becomes negative. So we intuitively 
guess that the shell has eigenmodes with negative lu 2 . 
Our numerical calculations show, however, that there is 
no such eigenmode. Since the shell we are considering 
couples to the bulk fields, it is considered that the lost 
of the energy by moving the shell is compensated by the 
increase of the bulk fields energy. 

Fig.2 
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FIG. 2. The typical configuration of the perturbed func- 
tions with the parameter g = 0.1, N = 100, V = 3000. The 
solid and the dotted lines are 8Z2 and 8Z§, respectively. The 
enhancon radius of these parameters is r e — 10.8123. The 
radius of the shell is ro = 20.0 and the eigenvalue is uj 2 = 1.0. 



Fig.3 
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FIG. 3. The typical configuration of the perturbed func- 
tions with the parameter g - 0.1, N - 100, V = 3000. 
The radius of the shell is ro = 10.8 and the eigenvalue is 
uj 2 — 1.0. The solid and the dotted lines are 8Z2 and 5Z$, 
respectively. Since the shell locates near the resonance radius 
Treso = 10.8071 the amplitudes of the perturbed functions 
become very large at the shell. 



VI. CONCLUSION AND DISCUSSIONS 

We have investigated the stability of the shell of 
wrapped D6-branes on K3 from the point of view of su- 
pergravity. As shown in Sec. [Il| if a probe D6-brane 
(wrapped on K3) is kicked inward slightly at the outside 
of the enhangon radius r — r e , the velocity reaches the 
speed of light at the radius. This behavior seems to be 
in contradiction to the behavior of the enhangon mecha- 
nism, which states that the D6-branes shell sits around 
the enhangon radius. So, what mechanism stabilizes the 
non-probe D6-brane? 

In the probe case, the background spacetime is fixed 
as a supersymmetric D6-D2 brane supergravity solution. 
Therefore, there is no obstacle to stop the motion of the 
probe since no force acts on it. In the non- probe case, 
the background spacetime gets excited by the motion of 
the shell and becomes non-supersymmetric. Then, new 
forces act on the D6 branes shell, as seen in Eq. ( 5.43[ ). 
In other words, the kinematic energy of the shell could 
be lost through the waves such as gravitational and the 
dilaton waves. This is one of the reasons why the non- 
probe D6-branes shell is stable. 

Quite interestingly, it is found numerically that there is 
no eigenmode with negative ui 2 for any radius of the shell. 
One may doubt it when the shell has negative tension, 
because its energy falls as it moves. This naive picture is 
not applicable to our case because the bulk energy is also 
changed. For example, let us consider the shell located 
inside of the enhangon radius. Because of the negative 
tension, the energy E falls as E oc — (SR) 2 when it moves 
slowly. On the other hand, the bulk energy of the dilaton 
field and the R-R fields may increase as the shell moves. 
So, we cannot discuss the stability of the whole system 
in terms of E only. Indeed, the D6-D2 system under con- 
sideration provides a notable example for demonstrating 
that D-branes shell with negative tension is stable in the 
sense that there is no eigenmode with negative to 2 . 

In Sec. V, we found that there is a resonant frequency 
mode whenever the shell is located at ro < r e . Because 
every time-dependent function /(£) contains this mode as 
a Fourier mode, the shell would oscillate with huge am- 
plitude. In this sense the shell is "unstable" gcncrically 
whenever the shell is located inside the enhangon radius. 
This fact suggests that the shell located inside of the en- 
hangon radius is pushed back to the enhangon radius, 
where it is stable. Therefore, we may say that the en- 
hangon radius is naturally selected as a suitable position 
of the wrapped D6-branes shell even from the viewpoint 



of the "classical" supergravity with the D-branes. This 
coincides with the enhangon picture stated in the intro- 
duction. 

It is interesting to investigate the stability of other 
models such as brane world scenario pa. There is a 
common belief that the brane with negative tension is 
unstable against the perturbation. By analogy with the 
present analysis, we may expect that it can be stabilized 
if one considers the perturbation of the whole system in- 
cluding the self gravity of the brane. It is also interesting 
to investigate the non-commutative effect of N D6-branes 
or non-linear effect of the D6-D2 branes system. After 
such steady works, we could see new physics of gravity 
combined with string theory. 
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